ARB I TRARY FUNCT ION GENE RAT ING CI RCUIT USING AS IMPLE OPERAT I ONAL 
ELEMENT AND AN ENCRYPTION METHOD THEREFOR 



5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to an apparatus for 
simulating complex systems at a high rate of speed and also to 
10 a method of modulating signals for encryption. 

2. Description of the Related Art 

Numerical models for describing quantitative changes are 
a type of arbitrary function generator that produces specific 

15 functions based on the parameters provided. The models can be 
simulated at high speeds with an analog arithmetic circuit and, 
in fact/ were simulated with analog computers in the past. 
However, these analog arithmetic circuits are rarely used anymore 
because the circuits have to be rewired for each problem and 

20 enormous circuits are required when the problems are of a large 
scale . Subsequently, wire connections in the analog arithmetic 
circuits were modified electrically through software, but this 
analog calculation method is no longer used because integrating 
complex circuits is extremely difficult. 

25 Inorder to simulate aphenomenonhavingmultiple variables 

at a high rate of speed, it is advantageous to generate arbitrary 
functions with an analog arithmetic circuit. For this method 
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to be practical, however, it is necessary to programmably 
generate functions of different types without needing to 
physically change the wiring of the arithmetic circuit . Further, 
the circuit must have a simple construction such as RAM that 

5 facilitates high integration by repeating simple mask patterns . 
It is also necessary to minimize the number of components needed 
to construct the circuit. 

A generalized Lotka-Volterra equation can be used in plant 
growth models or for the dynamics of biological communities and 

10 is represented by a differential equation such as the following. 



15 elements of a system; jiij is a constant of interaction between 
the elements; and ri is an inherent constant for each element 
representing the growth rate of the population i. When n is 
sufficiently large, this equation can approximate a positive 
arbitrary continuous function of an arbitrary precision by 

20 varying n and Uij, it is possible to generate patterns of x ± that 
change in various ways . Accordingly, in addition to modeling 
phenomena, this calculation has a wide range of applications 
as an arbitrary waveform generator. The present invention 
focuses on the symmetry of this equation to perform high-speed 

25 emulation by performing the calculation in hardware or using 
parallel processing with multiple digital computers. 



Equation 1 




Here, x± is the size of a population i representing the 
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SUMMARY OF THE INVENTION 



In view of the foregoing, it is an object of the present 
invention to provide a simple calculation element f and 
5 particularly an arbitrary waveform generating circuit or random 
function generating circuit for solving the generalized 
Lotka-Volterra equation through an arithmetic circuit employing 
relatively simple circuit elements. It is another object of 
the present invention to provide anencryption/modulationmethod 
10 employing the arbitrary function generating circuit. 

These objects and others will be attained by an arbitrary 
function generating circuit employing a simple arithmetic 
circuit element for calculating the following generalized 
Lotka-Volterra equation (Equation 1) . 
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Equation 1 — 
dt 



Here, x ± is the size of a population i representing the elements 
of a system; \Xx^ is a constant of interaction between the elements; 

20 and ri is an intrinsic constant for each element representing 
the growth rate of the population i. The arbitrary function 
generating circuit includes a plurality n of modules, each having 
a plurality n of input terminals and one output terminal; and 
a plurality n of connecting wires including a first connecting 

25 wire for connecting the output terminal of the first module to 
the first input terminal of each module; a second connecting 
wire for connecting the output terminal of the second module 
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to the second input terminal of each module; and continuing in 
a like manner until the nth connecting wire for connecting the 
output terminal of the nth module to the nth input terminal of 
each module, 

5 Each module comprises a group of n variable resistors 

connected to the n input terminals; an output sum connecting 
wire connecting the output terminals of the variable resistors 
in order to total their output values; a multiplier having an 
amplifier for multiplying the sum of the output values by the 

10 output value from the corresponding module; and an integrator 
for integrating the output values. 

The generalized Lotka-Volterra equation is converted to 
Equation 7 below by setting the interaction coefficient Uij 
equivalent to the value of the variable resistors in each module * 

15 An arbitrary module in the group of modules is designated as 
the first module to output constant value 1. A fixed value ri 
is replaced by the value in the variable resistors input 
via the first input terminal of eachmodule by the first connecting 
wire . 
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Equation 7 l^***'^***** ('^vO 



The arithmetic circuit calculates the converted 
generalized Lotka-Volterra equation by iterating the following 
25 process in which the size Xi (i = 1, 2, n) of a population 
i representing the elements comprising the system for generating 
an arbitrary function is output by the output terminal of each 
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of n modules. This output is transferred to a group of input 
terminals in a group of modules via n input connection wires 
and totaled by a connection wire. 

According to another aspect of the present invention, the 
5 arbitrary function generating circuit further includes an 
interaction constant setting means for programmably changing 
each value in the group of variable resistors in each module* 
Alternatively, FET semiconductor element circuits are provided 
as the variable resistors in each module . The resistance values 
10 in the variable resistors are varied externally using control 
signals . 

According to another aspect of the present invention, an 
arbitrary function generating circuit includes a plurality (n) 
of modules each having one input/output signal terminal, and 

15 one signal bus connecting each of the input/output signal 
terminals. Each module comprises a frequency synthesizer; a 
first multiplier for multiplying the output from the frequency 
synthesizer with the value input from the signal bus via the 
input/output terminal; a low pass filter for removing the AC 

20 component from the value output by the first multiplier; a second 
multiplier into which the output from the low pass filter is 
input; an integrator for integrating the output from the second 
multiplier; a connecting circuit for inputting the value output 
from the integrator into the second multiplier, such that the 

25 second multiplier can multiply this value from the value output 
by the low pass filter; an oscillator; and a third multiplier 
connected to the input/output terminal for multiplying the value 
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output by the oscillator with the value output from the 
integrator . 

The value xj. (i = 1, 2, n) output by the integrator 
in each module is the number of elements constructing a system 

5 for generating an arbitrary waveform function and is multiplied 
with a carrier wave signal having an angular frequency o> L 
generated by the oscillator by the third multiplier. The 
multiplied value XiSincoit is output from each module to the signal 
bus; and an addition value e as shown in Equation 8 is generated 

10 by multiplexing the output from each module and re-input via 
the signal bus . 

Equation 8 e - x a sin <u t t + x 2 sin to 2 1 + • • • + x m sin (O a t 

15 The value output from the frequency synthesizer is set 

to a voltage Wi shown in Equation 9 generated by multiplexing 
the interaction coefficient . 

Equation 9 W t - /j^sinco l t+(J. s siaQ} 2 t+---+fJ. lU siaoJ lt t 

20 

The addition value e is multiplied by the voltage W x in 
the first multiplier and passed through the low pass filter to 
obtain the DC component shown in Equation 10. 

25 Equation 10 eW ( - fi y x l +fi v x 2 +--+fM ni x ll 

This value eWi is multiplied by the value of Xi output 
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by the integrator in the second multiplier and integrated in 
the integrator to obtain the signal for Xi. The signal Xi is 
multiplied in the third multiplier by the carrier wave generated 
by the oscillator to obtain the signal XiSincoit. The signal 
5 Xisincoit is output onto the signal bus via the input/output 
terminal; and the process is repeated. 

According to another aspect of the present invention, the 
arbitrary function generating circuit further includes a means 
for setting the voltage Wi with the frequency synthesizer that 
10 is provided either external to the module as an amplifying circuit 
for extracting signals from the internal oscillator and adding 
weight to the signal or internal to the module as an independent 
circuit . 

According to another aspect of the present invention, an 
15 arbitrary function generating circuit includes a frequency 
synthesizer that outputs a value W; a first multiplier that 
receives the output value W as one input value; a first low pass 
filter for cutting out a first frequency component from the value 
output from the first multiplier; a second multiplier that 
20 receives the value e 4 output from the first low pass filter as 
one input value; a second low pass filter that removes a second 
frequency component lower than the first frequency component 
from the value output from the second multiplier; an adder that 
receives the output value e 5 as one input value; a delay circuit 
25 that delays the output value ei and sets the output value ei 
as the second input value to the adder; a first frequency 
multiplier that multiplies the frequency of the output value 
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ei and outputs an output value e 2 as the second input value for 
the second multiplier; and a second frequency multiplier that 
multiplies the frequency of the output value ei and outputs an 
output value e 3 as the second input value for the first multiplier . 
5 The amplitude of the signal x x (i = 1, 2, n) in Equation 

7 is expressed by the AC signal Xisin(itoo) . Here, i = 1, 2, 
.../ n. and coo is a basic angular frequency. The signal ei output 
from the adder is added to the total signal x ± to generate a 
frequency multiplexed signal as shown in Equation 11. 

10 

Equation 11 e l - jr 1 sin(o> 0 )+x 2 sin(2£e? 0 )f+-"+^sin(/iaj 0 )r 

The frequency multiplexed signal ei is constantly 
maintained in a closed loop delayed by the delay circuit only 

15 a time sufficiently longer than the period of the base frequency 
o>q. The first frequency multiplier multiplies the frequency 
of the signal ei by a (n « a ) and outputs the signal e 2 . The 
second frequency multiplier multiplies the frequency of the 
signal ei by j3 ( a « j3 ) and outputs the signal e 3 . The frequency 

20 synthesizer generates at once a signal W that multiplexes the 
interaction coefficient Uij (i and j — 1, 2, n) and outputs 
the signal W. The first multiplier multiplies the signals e 3 
and the W, The first low pass filter removes the frequency 
component greater than a frequency near j3co 0 and outputs the 

25 following signal e 4 . The second multiplier multiplies the 
signals e 2 and e 4 . The second low pass filter removes the 
frequency component greater than a frequency near a co 0 and outputs 
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the following signal e 5 shown by Equation 18. 



Equation 18 
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The generated signal e 5 , which corresponds to the frequency 
multiplexed differential signal dei/dt on the right side of 
Equation 7, is added to the signal ei in the adder to perform 
the integration of Equation 7. 

10 According to another aspect of the present invention, an 

encryption method using the arbitrary function generating 
circuit described above sets the initial value ei I t =o of the 
frequency multiplexed signal ei as the code to be encrypted and 
the signal W output from the frequency synthesizer as the key 

15 code for encryption. A time change pattern of the signal e x 
or the signal eilt=Tat a time T is encrypted. 

According to another aspect of the present invention, the 
signal W output from the frequency synthesizer is used as the 
code to be encrypted and the initial value ei I t =o of the frequency 

20 multiplexed signal e x is used as the key code for encryption . 



BRIEF DESCRIPTION OF THE DRAWINGS 



In the drawings : 

25 Fig. 1 is a block diagram showing a module of the first 

embodiment used in a random function generating circuit of the 
present invention; 
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Fig. 2 is a block diagram showing the overall schematic 
of an arithmetic circuit of the first embodiment; 

Fig. 3 is a block diagram showing a module of the second 
embodiment used in a random function generating circuit of the 
5 present invention; 

Fig. 4 is a block diagram showing the overall schematic 
of an arithmetic circuit of the second embodiment; 

Fig. 5 is a block diagram showing the overall schematic 
of an arithmetic circuit of the third embodiment having only 
10 one module; 

Fig. 6 is an explanatory diagram showing a decryption 
apparatus using the circuit of the present invention; 

Fig. 7 is an explanatory diagram showing paths indicating 
changes in the value of Xi over time on a plane formed by the 
15 two module outputs; and 

Fig. 8 is an explanatory diagram showing a public key 
encryption method using the system of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

20 

An arbitrary function generating circuit according to 
preferred embodiments of the present invention will be described 
while referring to the accompanying drawings. A first 
embodiment will be described with reference to Figs. 1 and 2. 
25 Fig. 2 shows an overall schematic for an arithmetic circuit of 
a generalized Lotka-Volterra equation, while Fig. 1 is a block 
diagram of a module in that circuit. 



Here, M ± represents an arbitrary module . The arithmetic 
circuit is configured of n number of modules. Each module has 
n number of input terminals Pij and one output terminal Q x . 

The arithmetic circuit further comprises a first 
5 connecting wire for connecting the output terminal Qi of the 
first module Mi to the first input terminals Pn, P21/ P n i of 
each module; a second connecting wire for connecting the output 
terminal Q 2 of the second module M 2 to the second input terminals 
P12/ ^22f -f p n2 of each module; and continuing in this way to 

10 an nth connecting wire for connecting the output terminal Q n 
of the nth module M n to the nth input terminal Pi n , F 2 nr .../ Pnn 
of each module. 

Each module further comprises n input terminals P X j, n 
variable resistors Vij connected to the corresponding input 

15 terminals Pij, connecting wires that add up the outputs from 
each of the variable resistors Vij , an amplifier A for amplifying 
the these totaled values, a multiplier T for multiplying the 
totaled value with the module output value, and an integrator 
S for integrating this output. 

20 The generalized Lotka-Volterra equation can be computed 

by connecting arbitrary modules Mi of Fig. 1 in the configuration 
shown in Fig. 2. Here, we will assume that an interaction 
coefficient Uij is equivalent to the value of the variable 
resistors Vij in each module and that one special module that 

25 outputs a fixed value is provided such that r x is the input value . 
Accordingly, Equation 1 can be rewritten as follows. 
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Equation 2 £l . ^ + ^ + J< ^ j (i -1,2,-y) 

If the first module M x is a fixed output module that always 
outputs the value Xi = 1, then we have the following. 

Equation 3 ±L mX ^ + ^ + g j (z=2,3/*v0 

Such that, 



Equation 4 Xi = 1 

Here, ri provides a fixed value (biased value) in the 
parentheses on the right side of the equation. If the role of 
ri is assumed by the weighted input nn from the first module 
Mi (xi = 1) , which constantly outputs a fixed value of 1, then 
by assuming fin is equivalent to r if the following equation can 
be formed. 

Equation 5 *i|ttt*i + Jj A***/j (i-2A—,n) 

Based on this assumption, 
Equation 6 = r ± 

By returning the term jinxi in the parentheses to the S on the 
right side of Equation 5, we have the following equation. 
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Equation 7 _L , ^ g (i-VvO 

Therefore, by introducing a module with a fixed output, 
Equation 1 can be simplified to Equation 7. The right side of 

5 the equation can be calculated with an adding circuit for 
performing a sum calculation (configured by the variable 
resistors V^ [j = 1-n] and the amplifier A in Fig. 1) and a 
multiplying circuit (the multiplier T in Fig. 1) for multiplying 
the value x ± by the result of the previous calculation. The 

10 values for Xi, when i = 2, 3, n, are obtained by integrating 
this result using an integrating circuit (the integrator S in 
Fig. 1) . 

The time change pattern for values of x if when i = 2, 3, 
n, can be changed arbitrarily by the value set for the 
15 interaction coefficient u iD (the value set in the variable 
resistors Vij [j = 1-n] in Fig. 1) . Therefore, the circuit 
connecting n number of modules Mi in Fig. 1 can be used as an 
n-1 channel arbitrary function generator. The value of Uij of 
the Vij can be varied externally using an element such as a field 
20 effect transistor (FET) . 

This circuit can generate functions having considerable 
complexity with several to several tens of modules . As the number 
of elements, or n in Equation 7, becomes large, the number of 
connections increases on the order n 2 , making integration 
25 difficult. An additional complex circuit must be provided 
separately to control the value pi] of the v in externally. 

Next, a second embodiment of the present invention will 
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be described with reference to Figs. 3 and 4. Fig. 4 shows the 
overall schematic of an arithmetic circuit for the generalized 
Lotka-Volterra equation, while Fig. 3 is a block diagram of the 
circuit modules in that arithmetic circuit . Here, Ni represents 
5 an arbitrary module. The arithmetic circuit is configured of 
n number of modules . Each module has one input/output terminal 
Ri. The input/output terminal R± for each module is connected 
to one signal bus line. 

Each arbitrary module N ± includes a frequency synthesizer 
10 Fl; a first multiplier Tl for multiplying the output from the 
frequency synthesizer Fl with the value input from the signal 
bus line via the input/output terminal R±; a low pass filter 
L for removing the AC component from the value output from the 
first multiplier Tl; a secondmultiplier T2 into which the output 
15 from the low pass filter L is input; an integrator SI for 
integrating the output from the second multiplier T2; a 
connecting circuit for inputting the value output from the 
integrator SI into the secondmultiplier T2, enabling the second 
multiplier T2 to multiply this value from the value output from 
20 the low pass filter L; an oscillator G; and a third multiplier 
T3 connected to the input/output terminal- Ri for multiplying 
the value output from the oscillator G with the value output 
from the integrator SI. 

When n is large, the input/output of the input/output 
25 terminals R± in the arbitrary modules Ni shown in Fig. 3 are 
set to sine waves of differing frequency (the amplitude is set 
as the output value) , and signal transfer is performed using 
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frequency multiplexing. With this method, all signals can be 
transferred over one connection line (signal bus) , achieving 
interconnection with a single connecting wire, as shown in Fig, 
4 . In other words, by setting the output value of the arbitrary 
5 module Ni to XiSincoit. Signals from each arbitrary module N± 
are added to one connecting wire. Accordingly, the voltage of 
the connecting wire is calculated as follows* 

Equation 8 e =*x x sin a> t t +x 2 sin co 2 t + x n sin co H t 

10 

The interaction coefficient ji^ is similarly expressed 
by an AC signal. Accordingly, an interaction constant is 
supplied to the x± as the following voltage, 

15 Equation 9 ^ - p y sinea x t+ij. 7J $inca 2 £+'*-+fj.^smc0 /t t 

This is generated independently by the frequency 
synthesizer Fl provided internally in each arbitrary module Ni 
or externally. The frequency synthesizer Fl canbe an amplifying 

20 circuit that adds weight to the signal output from the oscillator 
G disposed inside each module or can be provided in the module 
shown in Fig. 3 as a completely independent circuit in order 
to eliminate all problems of phase deviation. 

With the configuration of the arbitrary module Ni shown 

25 in Fig. 3, the first multiplier Tl generates a signal by 
multiplying e by Wi. This multiplied signal contains a DC 
component (a component that changes slower than the carrier wave 

15 
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sincoj.t) and an AC component (the frequency being the sum and 
product of each carrier wave frequency) . The DC component of 
this signal is obtained when the signal passes through the low 
pass filter L, as follows. 

5 

Equation 10 eW t - it u x x +fJL u x z + -~+iA td x n 

The sum on the right side of Equation 10 is obtained in 
real-time. The second multiplier T2 multiplies the value of 

10 Xi with this signal and outputs the result, which is the value 
on the right side of Equation 7. The integrator SI integrates 
this signal to obtain the signal for xi. Next, the third 
multiplier T3 multiplies the signal xi with a carrier wave sin 
coit signal generated by the oscillator G to generate an x±sin 

15 ojjt signal. This resulting signal is simultaneously input into 
the first multiplier Tl and output externally via the 
input/output terminal Ri. 

In this way, a calculation equivalent to the generalized 
Lotka-Volterra equation can be performed by connecting multiple 

20 modules N ± by a single signal wire, as shown in Fig. 4. 

A method for decreasing the number of connecting wires 
using different frequencies is a basic technology in analog 
multiplex communications. Reports for the hardware of neural 
networks (learning threshold elements) has already been 

25 published (Hirokazu Yokoi and Masao Saito, 1986; New Learning 
Elements-Foulethret; IEICE Transactions, Vol. J69-A, No. 6, 
pp. 1173-1175) . 

16 
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The second circuit type in the second embodiment of the 
present invention is novel in that it greatly reduces through 
signal f requencymultiplexing the number of wires between modules 
necessary for calculating the generalized Lotka-Volterra 
5 equation in a circuit. While the number of wires in the circuit 
of the first embodiment increases on the order of n 2 , that in 
the circuit of the second embodiment increases on the order of 
n . 

In the circuit shown in Fig, 4, n number of modules is 

10 sufficient to calculate a generalized Lotka-Volterra equation 
having n elements . Since a mo re complex function can be expressed 
with a larger number n of modules, it is desirable to reduce 
as much as possible the number of overall components and wires 
required for mounting. According to the next method, a 

15 generalized Lotka-Volterra equation of an arbitrary scale (an 
arbitrary n) can be calculated with only one module. 

Fig . 5 shows the third embodiment of the present invention . 
Fig. 5 is a block diagram showing the overall schematic of the 
arithmetic circuit for the generalized Lotka-Volterra equation 

20 achieved in one module. 

The module comprises a frequency synthesizer F2 that 
outputs a value W; a first multiplier T4 that receives the output 
value W as one input value; a first low pass filter LI for cutting 
out the first frequency component from the value output from 

25 the first multiplier T4; a second multiplier T5 that receives 
the value e4 output from the first low pass filter LI as one 
input value; a second low pass filter L2 that removes a second 
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frequency component lower than the first frequency component 
from the value output from the second multiplier T5; an adder 
B that receives the output value e 5 as one input value; a delay 
circuit V that delays the output value ei and sets the output 
5 value ei as the second input value to the adder B; a first frequency 
multiplier HI that multiplies the frequency of the output value 
ei and outputs an output value ei as the second input value for 
the second multiplier T5; and a second frequency multiplier H2 
that multiplies the frequency of the output value ei and outputs 
10 an output value e3 as the second input value for the first 
multiplier T4. 

By setting a carrier frequency according to a method 
described later andperf orming the calculation using the circuit 
shown in Fig. 5, it is possible to achieve the same effects as 
15 the circuit in the second embodiment with fewer components. 

The amplitude of the signal x ± , where i = 1, 2, n, in 
Equation 7 is expressed by the AC signal xisin(icoo) for i = 1, 
2, n. Here, o>o is the basic angular frequency and the value 
is arbitrary (the value is determined according to the objective 
20 and application and the frequency characteristics of the 
components during actual manufacturing) ♦ In other words, x x 
is expressed by the strength of the signal having an angular 
frequency of an integral multiple ofcoo. 

All signals can be added and transferred over one signal 
25 line, since each AC signal is linearly independent from the others 
(that is, signal transfer using frequency multiplexing) . A 
frequency multiplexed signal obtained by adding all frequencies 
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is shown in the following equation. 

Equation 11 e, -x 1 sin(a> 0 )+x 2 sin(2£t? 0 )t+—+x M sin(neo Q )t 

5 The frequency multiplexed signal ei is set as the input value 
of the delay circuit V in Fig. 5. 

The ei is delayed only a time sufficiently longer than 
the period of the base frequency co 0 (2tc/o>o) by the delay circuit 
Vand is again input into the delay circuit V after passing through 
10 the adder B, In other words, the ei is output from the delay 
circuit V and returns to the delay circuit V, and is therefore 
constantly maintained in this closed loop. 

The signal added to e x in the adder B is set to the frequency 
multiplexed differential signal dei/dt shown in the following 
15 equation. This differential signal is multiplexed in the same 
way as the differential value dxi/dt defined on the left side 
of Equation 7 , where i = 1, 2, n. 

de x dx+ . / \ dx. . u \ dx / \ 
Equation 12 ^ " ^HiK /+^m(2^ 0 >+ ..M-^in(*a> 0 > 

20 

Hence, ei and dei/dt are added by the adder B. 
Equation 13 

*i + ~- - ^ +^jsm(^ 0 >+^ 2 +^jsin(2oi 0 >+--^x I1 + ^jsin(/ia> 0 > 

Since this addition equation is exactly the integral 
calculationperf ormed according to the Euler method, the integral 
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calculation of Equation 7 is achieved by a closed loop circuit 
formed by the adder B and delay circuit V. 

Amajor issue is how to construct the frequency multiplexed 
differential signal dei/dt with a simple circuit. The present 

5 invention solves this problem by multiplying two signals having 
differing frequencies and twice making use of the phenomenon 
in which a beat signal corresponding to the differential 
frequency between the two is generated, as described below. 
The first frequency multiplier HI multiplies the frequency 

10 of the signal ei by a (n « a) and outputs the signal e 2 . 

Equation 14 e 2 - x l sm(aca Q )t+x 2 sin(2acQ Q )t+ ~-+x m Wi{naa> 0 ) 

Similarly, the second frequency multiplier H2 multiplies 
15 the frequency of the signal ei by j3 ( a « j3 ) and outputs the 
signal e3 . 

Equation 15 * 3 - x x ^{fi<o Q )t+x 2 m{2pai Q )t^-^x 9 wiflpa} Q ^ 

20 The frequency synthesizer F2 generates at once a signal 

W that multiplexes the interaction coefficient Uij, such that 
i and j = 1, 2, n, as follows, and outputs the signal W. 



Equation 16 

+ fi 2X sin(pco 0 + 2aco Q + 2co 9 )r+ sin(20a? u + 2aa> u + 2co 0 )+ * • sin(n£<u„ + 2a^ 0 + 2a> Q ) 

M 

+ H ni sin(pca Q + naco Q + /ia>„ )+;u fl2 sin(2£cE) 0 + naa; 0 + nm q )+ - -+ju w sin(/i£a>(, + /iaa; a + nco Q )r 
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With the construction described above, the first 
multiplier T4 first multiplies the e 3 and the W generated by 
the frequency synthesizer F2 to generate a signal comprising 
the sum frequency and difference frequency of an AC signal that 
5 includes e 3 and W. The first low pass filter LI removes the 
frequency component greater than a frequency near j3 o> 0 and outputs 
the following signal e 4 . 

Equation 17 

10 * 4 - | g jsin{a*> 0 + <o Q >+ | Jj ^ jsin(2aa;,+ 2a*. >+ • • m||| ^m{n*a>,+ no>« > 

Next, the second multiplier T5 multiplies e 2 and e 4 . The 
second low pass filter L2 removes the frequency component greater 
than a frequency near a co 0 and outputs the following signal es* 
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e 5 = x, 



This equation expresses the right side of Equation 2 that 
20 is multiplexed. In other words, this equation corresponds to 
the multiplexed differential signal dei/dt. 

The signal e$ (dei/dt) is added to the signal ei that has 
passed through the delay circuit V by the adder B, thereby 
integrating Equation 2 by the method described above. 
25 This method theoretically enables Equation 7 to be 

calculated for an arbitrary n using only one circuit shown in 
Fig. 5. Due to the difference in initial values of lil t =o, and 
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the signal W in Equation 16 in the circuits of Figs. 5, l t 
demonstrates various differing time change patterns. 
Accordingly, by interpreting li I t-o as the value to be encoded 
(modulated) and W as the key code for encryption, then the time 
5 change pattern of 1 ± or the signal 1± i t =T in time T can be used 
as the signal to be encoded (modulated) . It is also possible 
to assume the opposite, that is, that W is the value to be encoded 
(modulated) and li I t =o is the key code for encryption. 

As described above, it is possible to encode or modulate 
10 signals using the arithmetic circuits of Figs. 5. It is also 
possible to decode or demodulate the encoded signals to the 
original signal according to a key code using the same arithmetic 
circuits • 

Next, the signal decryption method will be described in 
15 greater detail. 

Fig. 6 shows an example configuration of a decryption 
apparatus 60 used for decoding encoded signals . xi canbe various 
functions depending on the value set for the interaction 
coefficient^, such that i and j = 1, 2, 3, n, and the initial 
20 value of x iA such that i = 1, 2, 3, n. Encryption uses this 
property as described below. 

The output values of n units are variables in which n-1 
variables change according to time, while one unit outputs a 
fixed value. For purposes of simplifying the description, we 
25 will choose two variables x k and x x and consider a path traveled 
by a point p at the coordinates x k and x x in the plane passing 
through these two variables (see Fig. 7) . According to the work 
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of the arbitrary function generation of the circuit, the path 
of the point p moves in an extremely large variety of patterns 
from an initial value A Point. Here we will assume that the 
interaction coefficient ]i± lf such that i and j = 1, 2, 3, 
n , converges at a point 1 when starting from the initial value 
A Point after describing a spiral path and converges at a point 

0 when starting from initial points B and C. We will further 
assume that the points of convergence are only these two points 

1 and 0. Accordingly, the transmitter can safely transmit 
one-bit data indicating 1 or 0. 

Hence, the transmitter can send a signal corresponding 
to the point A when wishing to send the data "1" safely and a 
signal corresponding to the point B when wishing to send the 
data "0." The receiver can accurately learn the data at the 
convergence point by applying the previously prepared 
interaction coefficient Uij to the decryption apparatus 60 . For 
example, when the transmitter sends data for the point A, the 
receiver inputs data for the point A and data for the interaction 
coefficient Uij into the decryption apparatus 60 to accurately 
learn that the convergence point is 1. 

Since the combinations of the interaction coefficient 
are nearly infinite, it is nearly impossible for someone unaware 
of the true value of the interaction coefficient to guess the 
interaction coefficient ]i±^ from the transmitted data. 
Accordingly, it would be extremely difficult to decode data that 
has been encoded by the decryption apparatus 60. 

Next, the method of encoding the signals with an encryption 
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apparatus will be described in more detail. 

When not using a public key (interaction coefficient y-ij ) 
The process below is followed in order to encode data 

equivalent to the point 1. 

(1) An initial value is generated arbitrarily. If the value 
converges to the point 1 in the encryption apparatus , that initial 
value is used as the encoded value. If, however, the initial 
value converges at point 0, that initial value is discarded and 
another initial value is generated arbitrarily. 

(2) Step (1) is repeated until an initial value that converges 
at point 1 is obtained. The procedure for encoding a value 
corresponding to the point 0 can be performed according to the 
same method. 

When using a partially public key 

When transferring encoded credit card numbers for mail 
orders and the like, it is necessary to transfer a key for the 
encryption by mail . This is highly inconvenient. Therefore, 
a public key encryption method is necessary. The conventional 
method using prime numbers utilizes the fact that an extremely 
large method would be required to factor out prime factors from 
large numbers. Accordingly, the method is flawed in that 
security drops abruptly with improvements in computer 
performance and the development of efficient prime factorization 
algorithms. The method described below encodes while 
performing simulation using analog operations, and therefore 
has extremely tight security, since the amount of calculations 
required for a third party to decode the data is nearly infinite. 
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The procedure for performing encoded data transfers by a public 
key encryption using the encryption apparatus is ciescribedbelow . 
Preparations 

A group of n units forming the present apparatus is divided 
into an A group and a B group. The A group is a group of units 
for the encryption apparatus possessed by the encoder and 
transmitter. The B group is a group of units for the decryption 
apparatus possessed by the decoder. Most of the interaction 
coefficient between the A and B groups of units are 0, while 
only a specific interaction coefficient is a value other than 
0. Only some of the units in the B group are connected by an 
interaction coefficient other than 0 to the A group and affect 
the B group. This group is called the C group. Similarly, a 
portion of units in the A group are connected to the B group 
by an interaction coefficient other than 0 and affect the A group. 
This group is called the D group. 

For simplicity, let us assume that the interaction 
coefficients are set to have only two convergence points, such 
that the output values x k/ xi, x a , and x b from four units have 
two convergence points in a four-dimensional area. The units 
k and 1 outputting the variables x k and xi belong to the A group 
(however, the output values do not belong in the public C group) 
and the units a and b outputting the variables x a and x b belong 
to the B group (however, their output values do not belong in 
the public D group) . The convergence point in the plane of 
variables x k and x x has a one-on-one correspondence with the 
convergence point in the plane of variables x a and x b . The 
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interaction coefficient and growth rate (fixed parameter r L ) 
are preset to enable estimation of the convergence point in the 
plane of x k and Xi from the convergence point from the plane 
of x a and x b . 

5 While the transmitter only possesses data concerning the 

A group, the transmitter can dynamically extract data regarding 
the D group from the receiver to perform calculations for the 
A group units. In other words, the transmitter can simulate 
(calculate) movement of the coordinates (x k , xj from data on 

10 the A group and D group. Similarly, the receiver can simulate 
movement of the coordinates (x a , x b ) from data concerning the 
B and C groups. 

Public key data 

To encode and transmit data, the transmitter acquires data 
15 of the interaction coefficient Uij and fixed parameter ri for 
the A group made public by the receiver and acquires of data 
of the output value x ± for the C group (i being the unit number 
of the C group) as unencrypted data. Using this data, the 
transmitter performs a simulation of the A group. 
20 Encryption accompanying normal interactive 

communications 

Since output values from the C group vary during the 
simulation, the transmitter simulates the Agroup while receiving 
unencrypted output values for the C group in real-time, in order 
25 to simulate the Agroup . At the same time, the receiver simulates 
the B group while receiving unencrypted data for the D group 
from the transmitter in real-time . In other words, the receiver 
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and transmitter exchange unencrypted data on the D and C groups 
in real-time, while the transmitter simulates the A group and 
the receiver simulates the B group independently. It is obvious 
that a third party cannot participate in this interactive 
communication . 

Data that the transmitter wishes to send encrypted and 
the initial values of the A group are not known to the receiver 
or a third party. The receiver can estimate data concerning 
the convergence value of the A group from position data for the 
convergence point in the B group simulation. However, a third 
party cannot simulate either the A or B groups because the third 
party does not possess the initial values for the A group or 
data for the B group. Accordingly, a third party cannot decode 
the encrypted data. 

The arbitrary function generating circuit according to 
the present invention constructed of a simple calculation element 
and the encryption/modulation method using this circuit 
demonstrates the following effects. 

(1) The present invention is a high-speed processing device for 
performing simulations. That is, the present invention can be 
used to perform high-speed simulations of a generalized 
Lotka-Volterra equation, which is employed as a complex and 
dynamicbiological model for simulating suchprocesses as changes 
in the number of species in a biosystem and the growth of plants . 
The arbitrary m function generating function of the present 
invention can be applied to a general simulator for such complex 
systems as market conditions and economic fluctuations and has 
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a wide range of applications as a high-speed simulator for complex 
systems . 

(2) The present invention can be used as an apparatus for 
generating complex waveforms. Accordingly, extremely complex 

5 waveforms can be generated with a relatively simple circuit 
simply by changing a small portion of the circuit parameters 
slightly, it is possible to create various waveforms. Since 
periodic or aperiodic waveforms can be generated according to 
the parameters, the invention can be used as a sound generator 
10 for electronic instruments or game machines. Further, since 
a slight change in parameters can generate irregular waveforms, 
the present invention can also apply to warning sirens and the 
like for preventing the threat of harmful pests. 

(3) The present invention can apply to an encryption apparatus. 
15 A generalized Lotka-Volterra equation has a plurality of 

equilibrium points according to the parameters and converges 
at different equilibrium points depending on the initial value . 
The set of initial values is infinite while the number of 
equilibrium points is finite. Therefore, data for secure 

20 transfer is encrypted as position data of a equilibrium point . 
When the transmitter sends the data, one initial value from the 
set of initial values that converge at this equilibrium point 
is arbitrarily selected. Using a circuit with the same parameter 
values, the receiver can find the convergence point. However, 

25 it is nearly impossible for a third party to decode the message 
without knowledge of these parameters. Accordingly, the 
present invention can be applied to an encryption apparatus for 
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high-speed and extremely secure data transfer. 

(4) The present invention can apply to a pattern recognition 
apparatus. Using a circuit having the same plurality of 
equilibrium points as the encryption apparatus described above, 

5 an initial value vector is set to the input pattern and a converging 
point is set to the recognition result. By establishing 
appropriate parameters in advance, the present invention can 
therefore be applied to a high-speed pattern recognition 
apparatus. Existing nonlinear optimal algorithms can be used 
10 to determine appropriate parameters . 

(5) The present invention can be used for encryption and 
decryption signals in spread spectrum communications. A 
circuit having a plurality of equilibrium points describes a 
certain path before arriving at the equilibrium point. This 

15 type of nonlinear system possesses robustness in achieving the 
same equilibrium point, even though a certain level of noise 
is generally introduced into the signal. The larger value n, 
the more robust the system is to noise . Accordingly, the circuit 
types 2 and 3 described in embodiments 2 and 3, respectively, 

20 having an n as large as possible are most appropriate. 

The output from the third multiplier of the type 2 and 
the output from the second frequency multiplier of the third 
circuit type have an extremely broad spectrum from the low 
frequency to high frequency range. After modulation, this 

25 signal can be transmitted in a transmission wave having an 
extremely broad band. The demodulated signal is input into a 
circuit on the receiving end having the same circuit constant 
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and the same components as those on the transmitting end (output 
from the third multiplier in the type 2 and output from the second 
frequency multiplier in the type 3) . The receiver can therefore 
obtain correct .data, providing the noise introduced into the 
5 signal is not large, since the signal describes the same path 
as that for the transmitter and reaches the correct equilibrium 
point. Even if the demodulated signal contains a large amount 
of noise, converging on atractor point will increase the 
percentage of the signal reaching the same equilibrium point 
10 as the transmitters. This percentage is greater with a larger 
n. 
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